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Abstract 



A method is presented for making finite Fourier mode truncations of the 
Rayleigh-Benard convection system that preserve invariants of the full partial 
differential equations in the dissipationless limit. These truncations are shown 
to have no unbounded solutions and provide a description of the thermal flux 
that has the correct limiting behavior in a steady-state. A particular low- 
order truncation (containing 7 modes) is selected and compared with the 6 
mode truncation of Howard and Krishnamurti,0 which does not conserve the 
total energy in the dissipationless limit. A numerical example is presented to 
compare the two truncations and study the effect of shear flow on thermal 
transport. 
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I. INTRODUCTION 



In a horizontal layer of fluid with fixed higher temperature on the bottom boundary and 
fixed lower temperature on the top boundary, cellular convective flow occurs for a certain 
range of Rayleigh number R and Prandtl number cr.l Such states of thermal convection are 
ubiquitous in nature, occuring in slightly modified form in the atmosphere, the ocean, the 
earth's mantle, and in the convection zone of the interior of stars. It was originally believed 
that flows in a finite container should scale as its vertical dimension. However, experiments 
have shown that thermal convection in a horizontal layer of fluid heated from below can 
show motions spanning the largest horizontal dimension of the container,! known as shear 
flows. Since the experimental setting precluded any externally imposed shear forces, it is 
concluded that these shear flows are driven by a Reynolds stress tensor with non-vanishing 
horizontal average. This behavior is also seen in numerical experiments.!^ 

In this paper we make truncations of the Rayleigh-Benard system using the Galerkin 
(or spectral) method.! The chosen basis is the standard Fourier one, because of its great 
simplicity and the fact that it is especially well-suited to the stress-free boundary conditions. 
One of the major points to be addressed is how to make a truncation that retains the 
truncated invariants of the full Partial Differential Equations (PDE's) in the dissipationless 
limit, and whether or not this has any effect on simulations in the dissipative case. We 
will show that energy- conserving truncations have several advantages over their non-energy- 
conserving counterpart. 

Section |lj is devoted to a presentation of the equations governing the system and the 
preserved quantities of the equations of motion in the dissipationless limit are displayed. 
In Section |T| we use the Galerkin method to expand the stream function and temperature 
field into a complete set of modes, deriving a set of coupled Ordinary Differential Equations 
(ODE's). The condition under which a mode truncation preserves the invariants of the 
dissipationless PDE's is obtained. Section [TV] presents a low-order model for shear flow 
generation, which is an extension of a model presented in Ref. [1]. The model is compared 
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with its predecessor and the advantages are explicitly shown. Finally, Section |V| is a summary 
of the main arguments of the paper. 



II. MODEL EQUATIONS 

The two-dimensional Rayleigh-Benard problem for an incompressible fluid is governed 
in the Boussinesq approximation by 



d\/ 2 tp 



+ 



i> , V 2 V>] = ^ + ^V 4 V, (la) 



Of 

dt Lr ' 1 dx 



^L + [i ,,T]= d 4 + ^T, (lb) 



where if) is the stream function, T is the deviation of the temperature from a linear conduction 
profile, v is the kinematic viscosity, and k is the thermal conductivity. All the quantities 
are dimensionless, the Prandtl number a is u/k, and the Rayleigh number R is 1/vk. The 
horizontal coordinate is x and the vertical one is y, with (x, y) £ [0, 2nL] x [0,n]. (The 
choice y £ [0, 7r], as opposed to y £ [0, 1], leads to a Rayleigh number smaller by a factor 
of 7r 4 .) The Poisson bracket used is defined as [A , B] = d x A d y B — d y A d x B, and the velocity 
field is given in terms of the stream function by v = z x Vip = (—d y ip, d x ip). We assume the 
fluid is periodic in x and has stress-free boundary conditions at the top and bottom walls: 

^ = V 2 ^ = d x ip = T = 0, for y = or vr. (2) 

In the dissipationless limit, v = k = 0, Eqs. ([!]) admit an infinite number of con- 
served quantities (this is a general feature of noncanonical infinite-dimensional Hamiltonian 
systems — for a discussion of the invariants of an analogous set of equations, see Ref. ^|). 
Here we shall concern ourselves with the total energy, E, given by 

E= l -{{Vi>Y)-(yT) = K + U, (3) 

where K is the kinetic energy, U is the potential energy, and the angle brackets denote the 
integral over the fluid domain. In the dissipationless limit, the time derivative of E is 
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E = -U 



dt Y dt I Y dt , 

= (ij([tP,V 2 ij]-d x T)) + (y 

= -(if>d x T) + ('<i> [T,y)) = 0, 



,T]-d x iP)) 



(4) 



showing that the total energy is conserved (we have set surface terms to zero in Eq. |4j). The 
total internal (thermal) energy of the fluid, (T), is also conserved. 



III. MODE EXPANSION 



In this Section we look at truncations of system (]!]) and their properties. In the first 
part we derive the equations of motion for modes of the truncations (Section [III A| ). In 
Section [III B| we examine the behavior of the invariants mentioned in Section |I| after they are 



truncated and show (in a manner similar to Ref. [7]) that they can be made to remain invariant 
by adding certain modes to the system. Such truncations will be called energy- conserving. 
Finally, in Section [111 C| we examine the important properties of these truncations, namely 



that the truncated system has no singular solutions and that the thermal flux is properly 
modeled. 



A. Derivation of the ODE's 

To turn the system of partial differential equations ([!]) into ordinary differential equations, 
we use the following normal mode expansions for the ip and T fields: 

ii>(x,V,t)= E i>mn(t)e« max+ny \ (5a) 

(m,n) eA^p 

T(x,y,t)= T mn (t)e^ max+n y\ (5b) 

(m,n)sAj< 

where a = 1/L is the inverse aspect ratio. The summations are over some sets A$ and At of 
modes (i.e., (m, n) pairs, where both m and n can be negative or zero). If both of these sets 
are infinite and contain all possible (m, n) pairs, then the equalities hold in @; otherwise, 
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the expansion is a truncation. Expansion (|5|) is more general than those used in Refs. HJ^J^ 
in two ways: first, it allows for a variable phase in the rolls (by allowing the Vw' s to be 
complex) and second, the expansion admits a non- vanishing shear flow part (the ipo n modes). 
The reality of the fields and the stress-free boundary conditions (|2|) lead to 

VW = Cm -n = -VV-n , (6a) 

T m n T_ m _ n T m _ n , (6b) 

so that if, say, ipu is in A^, then so are ^—1,1, and ip^i^i, and similarly for A?- 

Note that A^ and At need not contain the same modes. For convenience we de- 
fine p mn = a 2 m 2 + n 2 to be the eigenvalues of the operator —V 2 . If we insert Eqs. @ 
into the Boussinesq equations ([I]), we obtain the following set of coupled nonlinear ODE's: 



— I - - I - ■ (Xm T (7 \ 

Pmn 

a (m'n"-m"n') 1p m i n > %p m "n" , 

Pmn 



m +m — m 



^Tmn ^PmnTmn ic^Tf^^Pmn K^™) 

+ ^2 a (m'n"-m"n'} 1p m >n> T m ii n n . 

n'+n" = n 



B. Preservation of the Invariants 

The kinetic and potential energies of Eq. (|3|) have the expansions 

K = Pmn \4>mn\ 2 , (8a) 

m,n 

(-1) P 

t/ = iE — T »r (8b) 

We now ask whether the total energy E is still conserved in the dissipationless limit for a 
truncated system. Taking the time derivative of Eqs. (|S|) and using Eqs. (0) with v = k = 0, 
we obtain after some manipulation: 
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dt dV ' 

= ia J2(-^) P m^ mn T* m ^ p . (9) 

Let N = {max (n) | (m,n) e A^, U Ap}, i.e., the maximum vertical mode number included 
in the truncation. If we assume the sum over p runs from — N' to N', we can write Eq. (Rf) 



as 



d N N ' 

-E = iaY, E E (-l) p m^ mn T,; iri _ p . (10) 

ar m n=-Np=-N' 

Now replace p by s = n — p: 

j N N'+n 

/ = i«EE E (-ir™^r; sl (ii) 

aZ m n =-N s=-N'+n 

and note that the maximum lower bound for s is — N' + N when n = N, while the minimum 
upper bound is N' — N when n = —N. If N' = 2N, s G [— iV, iV] always (since for \s\ > N 
the mode is not included in the truncation and so is made to vanish), and we can use the 
symmetries given by Eq. (|6]) to show that dE/dt vanishes.@ Hence, we must have p running 
from — 2N to 2N, which from Eq. (|8b"| ) implies adding the modes T 0j _2at . . . T 0i 2tv to At- 
For the internal energy, (T), the expansion is 

" p odd V 

and its time derivative in the dissipationless limit is 

= ^T~E E m i» T m,»-r ( 13 ) 
al n rn,n p odd 

Comparing Eq. fll3|) with Eq. (§), we see that Eq. ( |13|) vanishes under the same condition 
as the total energy E% 



C. Properties of the Truncations 



To show that the truncated systems obtained in Section 
for all times t > 0, we consider the quantity 
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have bounded solutions 



Q = K + 2 £ \T mn \ 2 + ^ (TL (14) 

m,n>0 n>0 V (l ' 

where K is the kinetic energy defined previously. The quantity Q is non-negative and it 
includes all the modes in the truncation such that if any of them diverges, then Q diverges. 
Thus, if Q is bounded from above then the truncated system has no unbounded solutions. 
Taking the time derivative of Q, with the viscosity and thermal conductivity nonzero, we 
can write 

^-Q<- min{2z/, k}Q + 4«JV , (15) 

with N being the number of T | n | modes included in the truncation. For 
Q > 4:kNq/ mm{2u, k} } we have dQ/dt < 0, and so Q is bounded. 

We define the horizontally averaged vertical thermal flux as q(y) = qcv(y) + Qcd(y), 



where q C v(y) = v y T is the convective thermal flux and q c d(y) = k(1 — d y T) is the conductive 
thermal flux (the overbar denotes an average over x). For energy-conserving truncations one 
can write the expansion for q as0 

M = it) - 2 E^j t T L . (16) 

m>0 " L ai 

In a steady-state situation this reduces to the expected result q — (q), independent of y, 
showing that the energy cannot "pile up" in steady convection. For a general truncation 
(for example, the truncations in Refs. [3],|8],|9]) , one cannot write q in the form given by Eq. ( |I~6|) 
and the thermal flux has an unphysical y dependence in a steady-state. General truncations 
can also have unbounded solutions as is the case in Ref. [I] for large enough Rayleigh number. 

IV. LOW-ORDER TRUNCATIONS 

A popular truncation of the Boussinesq equations (in the spirit of the Lorenz modefl) 
is the 6-ODE model given by Howard and Krishnamurti^ and used by other authors.^ 
It includes the 6 independent modes ip l 0l1 ^n; ^12 > ^11 > ^12 > anc ^ ^02- The Howard and 
Krishnamurti truncation is the simplest one that allows for a nonzero shear flow (the "0oi 



mode). It has a vanishing (T). However, it is not energy- conserving: it lacks the Tq 4 mode. 
We will add this Tq 4 mode to the 6-ODE model to obtain what we will call the 7-ODE 
model. Figure |I| displays explicitly the energy-conserving property of the 7-ODE model in 
the dissipationless limit. Note that the nonconservation of energy for the 6-ODE model is 
not a small effect: Figure |I] shows that the energy is not conserved by a factor of 1.8 to as 
much as 2.8. 

Figure || shows the time-averaged quantity Nu x R/R c evaluated at the upper boundary 
plotted as a function of the Rayleigh number scaled by R c , the critical Rayleigh number 
where the fluid at rest becomes linearly unstable (R c = (1 + a 2 ) 3 jo? ~ 10.088 for a = 1.2). 
The Nusselt number Nu is the ratio of total heat transferred to the heat conducted when 
the fluid is at rest, so that Nu x R/R c is a dimensionless measure of the thermal flux.! 
The thermal flux in the 6-ODE model (solid line) is seen to grow rapidly after R/R c — 11, 
whereas the same quantity for the 7-ODE model has a slower growth (dashed line). The 
dotted line has a slope of 5.05, corresponding to the experimental result of Ref. |3] (the 
intercept is arbitrary) for o = 7 (the numerical results are also for a = 7). The agreement 
between the experimental results and that of the 7-ODE model is excellent. However, since 
the experiment was done with no-slip boundary conditions, caution should be taken in 
concluding the accuracy of the model (the value of R c used to scale the experimental curve 
is the no-slip one). The onset of shear flow (at R/R c — 8.1) is associated with a decrease in 
thermal flux for the 6-ODE model but an increase of that quantity for the 7-ODE model. 
For smaller values of a, the two models behave in a more similar fashion, as seen in Figure |3| 
where the same quantities are plotted for a = 1. 

Figure ^ shows the thermal flux q for the 6-ODE model evaluated at different values of y 
as a function of time (R/R c = 3.5, a = 1, a = 1.2). The thermal flux after the system 
settles in a steady-state is seen to depend on y, which is unphysical since this would lead to 
energy pile up (see Section [111 C| ). However, as shown in Figure ^ the thermal flux for the 
energy-conserving 7-ODE model is independent of y when the system reaches a steady-state. 
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V. CONCLUSIONS 



In this work we have developed a general method for generating energy-conserving 
Galerkin approximations of the PDE's that describe Rayleigh-Benard convection, using 
a generalization of the method of Treve and Manley! The truncations allow for shear flow 
(zonal flows independent of the horizontal coordinate, x) and variable phase of the rolls 
(breaking point symmetry with respect to the center of the rolls). These flows also have 
applications in tokamak plasmas, where it is thought that a shear flow in the edge layer 
is responsible for the so-called H-mode,ll! in which confinement is increased by a factor of 
two over the normal, or L-mode phase. Convection cells form as a result of the nonlinear 
development of the Rayleigh-Taylor instability in regions of unfavorable magnetic curvature. 
Such convection cell turbulence is widely observed in the edge of tokamak plasmas0 These 
vortices can lead to the generation of a shear flow in a manner analogous to the Rayleigh- 
Benard caseEHl§ and it is believed that this flow creates a barrier to particle transport, 
thereby improving confinement. The new truncations could help to provide a foundation for 



turbulence models of L-H transitions such as in Ref. [17. 

There are essentially three arguments for using energy-conserving approximations: First, 
the cascade of energy through the inertial range to the dissipation scale is modeled without 
extraneous terms in the energy equations. Thus, energy- conserving truncations are a good 
technique to reduce unphysical numerical dissipation or sources and possible instabilities. 
This makes the truncations more closely related to the full equations, since in the full 
PDE's the dissipation comes entirely from the linear terms and the nonlinear parts of the 
equations conserve energy. Note, however, that for the example given (the 7-ODE model) 
there are too few modes to speak of an energy cascade through an inertial range. The second 
property is the correct description of the thermal flux in the steady-state limit, even with 
dissipation. This means that the fact that a truncation does not conserve energy in the 
ideal limit significantly affects the way the energy flows in the dissipative regime, and so the 
energy-conserving truncations are relevant to the dissipative case. Finally, the boundedness 
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of solutions is a strong point since physically one does not expect divergent behavior and 
hence the energy-conserving truncations are more reliable. 

In the dissipationless limit, we have demonstrated that energy conservation is violated 
drastically for a typical truncation (the 6-ODE Howard and Krishnamurti model), whereas 
it is conserved to machine precision for the energy- conserving truncation (7-ODE model). 
In the presence of dissipation, the numerical results for the comparison of the two models 
for Prandtl number o = 7 clearly show that the energy- conserving truncation is much closer 
to experimental results. If the system reaches a steady-state, the 7-ODE model properly 
models the thermal flux, whereas the 6-ODE model exhibits an unphysical dependence on 
the vertical coordinate. Thus we believe that energy-conserving truncations represent the 
full system more accurately. Note, however, that the 7-ODE model is presented here as an 
illustration of the energy-conserving truncation technique, as severely truncated systems may 
yield numerical results far removed from the behavior of the full system.S'Uli For instance, 
in Ref. |]the authors concluded that there are not enough modes in the 6-ODE truncation to 
adequately model the bifurcations observed in a full simulation of the flow. Note also that 
Ref. [TP] gives a different 7-ODE extension of the 6-ODE model that conserves total vorticity 
in the dissipationless limit. That model could be made energy-conserving by adding to it 
a Tq 6 mode. 
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FIGURES 

FIG. 1. Total energy E in the dissipationless limit (y = k = 0) for the 6-ODE (solid line) 
and 7-ODE (dashed line) models with a = 1.2 and initial conditions (a) ipli = ipoi = ^12 = ^11 = 
all other modes begin at 10 -5 , (b) ip^ = T[ 1 = 1; all other modes begin at 10~ 8 . 

FIG. 2. Time-averaged Nu x R/R c evaluated at y = tt as a function of R/R c for the 6-ODE 
(solid line, triangles) and 7-ODE (dashed line, black dots) models. The dotted line has a slope 
of 5.05, corresponding to the experimental results for a = 7 in Ref. |3[ 

FIG. 3. Time-averaged Nu x R/R c evaluated at y = tt as a function of R/R c for the 6-ODE 
(solid line, triangles) and 7-ODE (dashed line, black dots) models. 

FIG. 4. Thermal flux q for the 6-ODE model evaluated at different values of y as a function of 
time. The final, steady-state thermal flux is seen to depend on y. 

FIG. 5. Thermal flux q for the 7-ODE model evaluated at the same values of y as in Figure § as 
a function of time. The final, steady-state thermal flux is independent of y as expected physically. 
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Thiffeault and Horton, Fig. 2 
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Thiffeault and Horton, Fig. 3 
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Thiffeault and Horton, Fig. 4 
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